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1. Introduction 

Chemical abundance data for the Galaxy shows a wide dispersion. We 
try to check if this dispersion may also be found with chemical evolution 
models as an effect of variations of an Initial Mass Function (IMF) which 
follows a Poisson's distribution. 

We will use two different methods: a) an analysis of the model internal 
error and b) a statistic study on the effects of the IMF discretization, 
based on montecarlo models. Both methods have been used in a chemical 
evolution model applied to the Solar Neighborhood. Preliminary results 
relative to O and N are shown. 

2. The Initial Mass Function 

We will try to justify the data dispersion by the possible variations in 
the model results which could be produced by the behavior of the IMF, 
whose properties are described by a probability distribution function: 

If we assume that N^t stars are observed in a mass range [mo, m^ax]) 
the mass m^ of the i-th star is a random variable whose probability 
distribution function is given by the stellar initial mass function: <I>(mj), 
whose integration along the whole mass range is normalized to 1. 

The number Nj of stars of mass m, is another random variable whose 
distribution function follows a Poissonian behavior, with the value of 
the IMF at that mass as the only parameter of the distribution: drii = 
dNi/Ntot = ^{mi)dmi, since the total number of stars is: Ntot = / dNi. 

A simple test to check the Poissonian nature of the distribution of Uj 
is simply the ratio of its variance to its average value as a function of 
mass. This ratio should be close to 1. Cerviho et al. (2001) illustrates 
in their figure 3 the results of a test with 1000 Montecarlo simulations 
of clusters with 10^ and 10^ stars with a Salpeter IMF slope. Despite 
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Hj is the ratio of a Poisson variable with a constant, it is Poissonian 
distributed within a 10%. 
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Figure 1. Dispersion of O - panel a)- and N -panel b)- abundance data of the 
Galaxy taken from different sources 

Once we know that nj can be approximated by a Poisson variable, it 
is possible to apply a proper statistical formalism into the code of our 
chemical evolution model in order to study the fluctuations which we 
may obtain from this stoscastic behavior for the IMF. 

3. The theoretical multiphase chemical 
evolution model 

Our modeled physical system is a cylindrical sub-galactic zone cen- 
tered on the Solar Neighborhood, extending about 1 kpc in the Galactic 
plane. The model used in this work is the code from Ferrini et al. (1992), 
improved with recent yields (Gavilan ct al. 2002). It computes a pro- 
duction matrix Qjj(m), where each (i,j) element represents the fraction 
of the star mass initially in the form of chemical species j, transformed 
and ejected as chemical species i, for 15 chemical elements. A matrix is 
computed for each stellar mass wcigthcd by the IMF. The nucleosynthe- 
sis production from SNI, SNII, and normal stars is taken into account. 
We let the model evolves for 13,2 Gyrs, with a 0.01 time step. 

The IMF used is taken from Ferrini et al.(1991): 

$(m) = 2.0865m-0-52xl0-P-07(logm)2+1.921ogm+0.73]5 

We will use the two described above methods: a) an analysis of the 
model internal error and b) a statistic study on the effects of the IMF 
discretization, based on montecarlo models. 



3.1 Error analysis 



3 




Figure 2. Variability of the IMF as a function of mass. 

Considering the Poissonian nature for the IMF, we assume an initial 
variance of cj^(<I>) = <I>. Then, we trail the error, following error propa- 
gation rules, within the model in every equation involved and related to 
our IMF. As a result we can obtain a measure of the variability of our 
model along the whole mass interval considered. 

We may observe in Fig. 2 that the variability increases strongly in 
the mass interval from 4 to 8 M©. This is exactly the interval of mass 
producing Fc and N. Therefore, in agreement with this result, our model 
must show a larger dispersion in these elements in comparison with other 
elements as O. 

3.2 Monte Carlo Simulations 

We use an uniformly distributed random numbers generator to gen- 
erate from it, and using the accumulated probability distribution, a suc- 
cession of numbers with our IMF as probability function. Our aim is the 
computation of a complete set of Qi^j Monte Carlo simulations, where 
the IMF is randomly different for each of them, with masses from 0.8 
up to 100 M0, divided in 800 intervals with linear interpolation within 
them. We also would like to study the effects of the discretization caused 
by rounding in the code. Our aim is to calculate ~ 500 — 1000 simu- 
lations, with ~ 10^ — 10^ stars. We will repeat our simulations for two 
different round down values: 10~^ and 10~^. The well calibrated model 
from Gavilan et al.(2002) will be our base of mean values. 
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4. Results 

The preliminary results have been obtained by the realization of 400 
montecarlo simulations for the galaxy disk zone, without having into ac- 
count the rounding in the code and with non calibrated models. Despite 
of these facts, a dispersion of ^ 10~^ already appears in our O and N 
abundances, as we may see in Fig.3. This fact seems to indicate that a 
good percentage of the observed dispersion could be achieved with the 
complete set of simulations. 
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Figure 3. Preliminary results of dispersion in O and N abundances. 



5. Conclusions 

With this kind of calculations we might estimate the variability and 
the errors that can be obtained by the models. Thus, we could include 
into the code the dispersion appearing as a consequence of the uncer- 
tainties in the IMF, and this way we could check if is possible to get 
the same dispersion as observed. It seems, from our preliminary results, 
that taking into account the sampling fluctuations from the IMF into 
the models can reproduce the observed dispersion. 
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